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Abstract. A thermodynamically consistent two-phase Stefan problem with 
temperature-dependent surface tension and with or without kinetic under- 
cooling is studied. It is shown that these problems generate local semiflows 
in well-defined state manifolds. If a solution does not exhibit singularities, it 
is proved that it exists globally in time and converges towards an equilibrium 
of the problem. In addition, stability and instability of equilibria is studied. 
In particular, it is shown that multiple spheres of the same radius are unstable 
if surface heat capacity is small; however, if kinetic undercooling is absent, 
they are stable if surface heat capacity is sufficiently large. 



1. Introduction 



In the recent paper |18| the authors studied Stefan problems with surface tension 
and with or without kinetic undercoohng which are consistent with thermodynam- 
ics, in the sense that the total energy is preserved and the total entropy is strictly 
increasing along nonconstant smooth solutions. 

1. To formulate this problem, let C M" be a bounded domain of class C^, 
n > 2. is occupied by a material that can undergo phase changes: at time t, 
phase i occupies the subdomain ili{t) of il, respectively, with i = 1,2. We assume 
that dfli{t) n dfl = 0; this means that no boundary contact can occur. The closed 
compact hypersurface T{t) := dili{t) C forms the interface between the phases. 

The problem consists in finding a family of closed compact hypersurfaces T{t) 
contained in 57 and an appropriately smooth function u : M-|_ x f2 ^ R such that 

' K{u)dtu — div{d{u)Wu) = 

d^u = 

N = o 

l^j{u)j + aU ^ j{u)V 

ld{u)d,u] - {l{u) - j{u)V)V 
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Here u{t) denotes the (absolute) temperature, iy{t) the outer normal field of ili{t), 
V{t) the normal velocity of r(t), n{t) = n{T{t)) = -divr(t)i'(t) the sum of the 
principal curvatures, and |w] — W2|r(t) ~ ''^ilr(t) the jump of a function v across 
T{t). Since u means absolute temperature we always assume that u > 0. 
Several quantities are derived from the free energies ipii^) as follows: 

• ei{u) := ipi{u) + urii{u) denotes the internal energy in phase i, 

• rjiiu) :— —ip'^iu) the entropy, 

• Ki(u) :~ ([{u) ~ —u'ip'l{u) the heat capacity, 

• l{u) := u|?/;'(u)] — —u\ri{u)\ the latent heat. 

Furthermore, di{u) > denotes the coefRcient of heat conduction in Fourier's law, 
7(u) > the coefficient of kinetic undercooling, and a > the coefficient of surface 
tension. In the sequel we drop the index i, as there is no danger of confusion; we 
just keep in mind that the coefficients in the bulk depend on the phases. 

The temperature is assumed to be continuous across the interface. However, 
the free energy and the conductivities depend on the respective phases, and hence 
the jumps ip{u) := |7/;(u)], |k(m)], |??(m)], |rf(u)] are in general non-zero at the 
interface. Throughout we require that the heat capacities Ki{u) and diffusivities 
di(u) are strictly positive over the whole temperature range u > 0, and that ip has 
exactly one zero Um > called the melting temperature. 

If we assume that the coefficient of surface tension a is constant, then this 
model is consistent with the laws of thermodynamics. In fact, the total energy of 
the system is given by 

E(u,r)= / e{u)dx+ / ads, (1.2) 
Jn\r Jr 

and by the transport and surface transport theorem we have for smooth solutions 
^E(u(i), r{t)) = - J^{ld{u)d,uj + le{u)jV + anV} ds 

= - J^{ld{u)d,uj - {l{u) - ^{u)V))V} ds - 0, 

and thus, energy is conserved. Also the total entropy $(w,r) defined by 

$(M,r) = [ f]{u)dx (1.3) 
Jn\r 

is nondecreasing along smooth solutions, as 

^J>(u(t),r{t)) = f ^d{u)\\/u\^ dx ^ I -{{d{u)d,u\+uli^{u)lV}ds 

= [ ^d{u)\\/u\^dx+ I --f{u)V^ ds > 0. 
Jq m Jr u 

2. In this paper we consider the physically important case where surface tension 
a — a{u) is a function of surface temperature u. Then, following [9^ and [T], 
the surface energy will be eriu) ds instead of J^ads, where er(w) denotes the 
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density of surface energy. In addition, one has to take into account the total surface 
entropy J^r]r{u) ds, as well as balance of surface energy. The latter means that 
the Stefan law has to be replaced by a dynamic equation on the moving interface 
r(t) of the form 

Kr(u)dt,nu - divr(dr(M)Vru) = ld{u)d^uj ~ {l{u) ~ -f{u)V + lr{u)n)V, 

where dt^n denotes the time derivative in normal direction, see (|1.7p . As in the 
bulk we define on the interface 

• er(u) := ct(u) + urir{u) denotes the surface internal energy, 

• r]r(u) := —a'{u) the surface entropy, 

• Kr(u) '.— Crl"") ^ —ua"{u) the surface heat capacity, 

• It{u) :— uct'{u) — —uririu) the surface latent heat. 

We also employ a Fourier law on the interface to describe surface heat conduction, 
i.e. we set qr ■= — (ir(u)Vru, which should be present as soon as the interface 
has heat capacity. Recall that u is assumed to be continuous across the interface, 
hence the surface temperature 

ur := M|p (1.4) 

is well-defined. 

Obviously, if a is constant then e-p ~ a, and rjr — kt — It — 0, hence this model 
reduces to (jl.ip . On the other hand, if a is linear in u we still have Kr — and 
then it makes sense to also set dr = 0, to obtain the modified Stefan law 

ld{u)d,uj = {l{u) - j{u)V + lr{u)n)V, 

which differs from the Stefan law in only by replacing l(u) by l{u) + It(u)'H. 
This is just a minor modification of (jl.ip . and its analysis remains essentially the 
same as in [18] . The only difference is that the stability condition for the equilibria, 
and in case 7 = also the well-posedness condition, changes. More precisely, the 
well-posedness condition changes from 7^ to A' 7^ where A(s) := (p{s)/a{s), 
and the stability condition modifies by replacing (p' /cr by A'. 

Therefore we concentrate here on the case where ky{u), dr{u) > 0, which means 
that a is strictly concave. It has been shown experimentally that positive surface 
heat capacity Kr is important in certain practical situations; see for recent work 
in this direction. Experimental evidence also show that a is strictly decreasing, 
hence admits exactly one zero Uc > 0; cr(u) is positive in (0, Uc) and negative for 
u > Uc- Physically, it is reasonable to assume Uc > Um- It turns out that the 
analysis of the problem with nonlinear surface tension is considerably different 
from the linear case. In the sequel we always assume that 

d,,^i,dT,a,j e C^{0,Uc), di, Ki,dr, Kr,cr > on (0, Uc), « = 1,2, (1.5) 

if not stated otherwise. Further, we let 7 = if there is no undercooling, or 
7 > on (0,Mc) if undercooling is present, and we restrict our attention to the 
temperature range u e (0, Uc)- 
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With these restrictions on the parameter functions, we consider the following 
problem: 






u 



m 



K,{u)dtu — div((i(u)Vu) 

M = 0, UT 

ifiiur) + cr(ur)^ 
Kr{ur)dt^nUr ~ divr(dr('"r)VrMr) = 

= ld{u)d^uj - {l{ur) + lr{ur)n - -f{ur)V)V on 

u{0) 

m 



n\T{t) 



Uq 

To- 



on dil 
on r{t) 
on r{t) 



m 



(1.6) 



Here f{u) ~ lip{u)l, and dt^nUr denotes the time derivative of ut in normal 
direction, defined by 



dt,nUT{t,p) -.^ —ur{t + T,x{t + T,p))\^^^, t>0, peT{t), 



(1.7) 



with {x{t + T,p) E R" : {t,p) G (—£,£) x T{t)} the flow induced by the normal 
vector field {Vi'). That is, [r x{t + T,p)] defines for each p S r{t) a flow line 
through p with 



dr 

and x{t,p) 



X{t + T,p) ^ {Vv){t + T,x{t + T,p)), X{t + T,p) eT{t + T), T €(-£,£), 

p. The existence of a unique trajectory 

{x(t + T,p) e R" : r e (-£,£)}, er(i), 

with the above properties is not completely obvious, see for instance |12) for a 
proof. 

Note that the (non-degenerate) equilibria for this problem are the same as 
those for ()l.ip : the temperature is constant, and the disperse phase fii consists 
of finitely many nonintersecting balls of the same radius. We shall prove that 
such an equilibrium is stable in the state manifold SAA defined below if VLi is 
connected and the stability condition introduced in the next section holds. Such 
an equilibrium will be a local maximum of the total entropy, as we found before in 
[18] for the case of constant surface tension. To the best of our knowledge, there 
is no mathematical work on thermodynamically consistent Stefan problems with 
surface tension depending on the temperature. 

3. The case where undercooling is present is the simpler one, as both equations 
on the interface are dynamic equations. In particular, the Gibbs-Thomson identity 

^{ut)V — a{ur)H = fiur) 

can be understood as a mean curvature flow for the evolution of the surface, 
modified by physics. 
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If there is no undercooling, it is convenient to eliminate the time derivative of 
ur from the energy balance on the interface. In fact, differentiating the Gibbs- 
Thomson law w.r.t. time t and, with A(s) = ip{s)/a{s), we obtain 

X'{ur)dt,nur + n'{r)V = on r(i), 

hence substitution into surface energy balance yields with 

Triur) :== wr(wr) - WlT), wr(ur) := A'(-ur)(^('"r) - Zr(ur)A(Mr))/Kr(wr), 

the relation 

Tr{ur)V = ^^^{diyr{dr{urWrur) + ld{u)d,uj}. (1.8) 

As V should be determined only by the state of the system and should not depend 
on time derivatives of other variables, this indicates that the problem without 
undercooling is not well-posed if the operator Tr{ur) is not invertible in L2{T), as 
V might not be well-defined. On the other hand if Ti^{uy) is invertible, then 

V = [Triur)]-' ^^^^\dWridr{ur)S/rur) + ld{u)d,u]} (1.9) 

uniquely determines the interfacial velocity V, gaining two derivatives in space, 
and showing that the right hand side of surface energy balance is of lower order. 
Note that 

c^r(s) = sa{s)[\'is)]yKr{s) > in (0, u,), (1.10) 
and a;r(s) = if and only if A'(s) = 0. Therefore the well-posedness condition 
becomes more complex compared to the case kt = 0. 

Going one step further, taking the surface gradient of the Gibbs-Thomson re- 
lation yields the identity 

Kriur)V - dr{ur)n{r) = Kr(Mr){/r(Mr) + Friu, ur)}, (1.11) 

as will be shown in Section 6. Here the function fr is the antiderivative of 
X{dr / Kr)' vanishing at s = and -Fr is nonlocal in space and of lower order. So 
also in the case where undercooling is absent we obtain a mean curvature flow, 
modified by physics. 

We would like to point out a phenomenon, in absence of kinetic undercooling, 
which is due to positive surface heat capacity Kr- If kt at an equilibrium is large 
enough, then such a steady state is stable, even if the interface is disconnected! 
Hence, this case seems to prevent the onset of Ostwald ripening. However, as we 
shall see in the next section, such equilibria cannot be maxima of the total entropy. 

The plan for this paper is as follows. In Section 2 we discuss some fundamental 
physical properties of the Stefan problem with variable surface tension. In partic- 
ular, it is shown that the negative total entropy is a strict Lyapunov functional 
for the problem, and we characterize and analyze the equiliria of the system. The 
direct mapping method based on the Hanzawa transform, first introduced in |8| , is 
discussed in Section 3. This way the problem is reduced to a quasilinear parabolic 
problem. In Section 4 we consider the full linearization of the problem at a given 
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equilibrium, and we prove that these are normaUy hyperbohc, generically. The 
last two sections deal with the analysis of the nonlinear problem with and without 
kinetic undercooling. The analysis is based on results for abstract quasilinear par- 
abolic problems, in particular on the generalized principle of linearized stability, 
see [lOl [17] . We refer here to [3l [131 [E] for information on maximal regularity in 
Lp- and weighted Lp- spaces, and to [3 HI [Tj [18] for more background information 
concerning the Stefan problem. 

2. Energy, Entropy and Equilibria 

(a) The total energy of the system (|1.6p is given by 

E(M,r)= / e{u)dx+ / er(ur)rfs, (2.1) 
Jn\r Jr 

and by the transport and surface transport theorem we have for smooth solutions 
^E(u,r) = J K{u)dtu dx - J le{u)]V ds + J {Kr{ur)dt^nUr - eriur^V} ds 

= J{-ld{u)d,ul - le{u)lV + divr(dr(wr)Vr«r) 

+ - (l{u) + lr{ur)n)V + -f{ur)V^ - er{ur)HV} ds 

= - j { IVX'^)] + (y{ur)n - -i{ut)V}V ds ^ Q 

by the Gibbs-Thomson law, and thus, energy is conserved. 

(b) The total entropy of the system, given by 

$(M,r)= / 7^{u)dx+ / r/r(wr)ds, (2.2) 
Jn\r Jr 

satisfies 

^$(u,r) = ^ rj'iu)dtudx + J {dt,nririur) - Mu)] + rjriur)n)V} ds 




{Kr{ur)dt,nUr + {l{u) + lr{ur)'H)V} ds 



+ I — {-Id(u)a,u]+divr(dr(ur)Vrur) + [d(u)9.w]+7(ur)l^'}ds 
Jr Ur 

= / ^d{u)\\/u\^ dx + I ■\{dr{ur)\'^rurf +urj{ur)V^}ds>0, 
Jn Jr Wp 

where we employed the transport theorem, the surface transport theorem and 
(|1.6|) . In particular, the negative total entropy is a Lyapunov functional for prob- 
lem pr^ . 
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(c) Even more, — $ is a strict Lyapunov functional in the sense that it is strictly 
decreasing along smooth solutions which are non-constant in time. Indeed, if at 
some time to > we have 

^$(u(io),r(io)) = o, 

then 

/ -\d(u)\Vu\'^ dx + [ \d(u)\Vurf ds + [ —-f(ur)V^ds = 0, 
Jn Jr -'r "r 

hence Vu(to) = in 51 and 'y{ur{to))V{tQ) = on r(to)- This implies u{tQ) = 
const = ur{to) in fl, and 'H{ta) — —l'ip{u{to))l/(T{ur{ta)) — const, provided we 
have 

|^(s)l=0 ^ a{s)>0. (2.3) 
Physically, this assumption is plausible, as it means that at melting temperature 
Um > (defined as the unique positive zero of the function (p{s) := |V'(s)I) the 
surface tension a{um) is positive. Since 57 is bounded, we may conclude that r(to) 
is a union of finitely many, say m, disjoint spheres of equal radius, i.e. {u{tQ) ,T {to)) 
is an equilibrium. Therefore, the limit sets of solutions in the state manifold defined 
below are contained in the {mn + 1)- dimensional manifold of equilibria 

£ = {(«*, [J Sr,{xi)) -.Q <u^ <Uc, |-0(u*)] = (?^ - l)o'(u*)/i?*, 

l<l<rn (2.4) 

BrM) C rj, BR,{xi)nBR,{xk)=^. k^l}, 
where Sr.^ {xi ) denotes the sphere with radius i?* and center xi . 

(d) Another interesting observation is the following. Consider the critical points 
of the functional $(u,ur,r) with constraint E(u,ur,r) = Eq, say on 

C/:={(u,ur,r): ueC(n\r), VeMU^iVt), MreC(r), u,Mr>0}, 

see below for the definition of Ai'H^{^^). So here we do not assume from the 
beginning that u is continuous across F, and Mr denotes surface temperature. 
Then by the method of Lagrange multipliers, there is /i G R such that at a critical 
point (it,, Ur*, r*) we have 

$'(u,, ur», r,) + /xE'(u,, Ur*, r,) = 0. (2.5) 

The derivatives of the functionals are given by 

($'(u,ur,r)|(u,ur,/^)) = {il'{u)\v)n + (77M«r)|«r)r - {Uu)} + 7^r{uT)'H{T)\h)r , 
and 

(E'(u,ur,r)|(M,«r,/j)) - (e'(u)|u)j, + (ei.(ur)|«r)r - ([e(")l + er(wr)H(r)|/i)r. 
Setting first vy = h = and varying v in (j2.5p we obtain 

77'(m*) + ^e'(u,) = in ri, 

varying ur yields 

77r(ur*) + /xep(ur*) = on F*, 
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and finally varying h we get 

lv{u,)j + Tjr{ur*)n{r) + Ai([eK)l + er(ur*)H(r,)) = on T,. 

The relations ?7(u) — -~^/{u) ande(M) — tp{u)~wjy (u) imply — — 
and this shows that — — is constant in fl, since k{u) — —uip"{u) > for 
all M > by assumption. Similarly on we obtain ur* — constant as 

well, provided Kr(wr) > 0, hence in particular = ur*- This further implies the 
Gibbs- Thomson relation [■(/'(u*)] + (t{u^)'H{T^) = 0. Since is constant we see 
that ■H(r*) is constant, by (|2.3p . Therefore T* is a sphere whenever connected, and 
a union of finitely many disjoint spheres of equal size otherwise. Thus the critical 
points of the entropy functional for prescribed energy are precisely the equilibria 
of problem (|1.6p . 

(e) Going further, suppose we have an equilibrium e, := (M*,ur*jr^) where 
the total entropy has a local maximum w.r.t. the constraint E = Eq constant. 
Then I?» := [$ + /iE]"(e*) is negative semi-definite on the kernel of E'(e,), where 
= — l/ii* is the fixed Lagrange multiplier found above. The kernel of E'(e) is 
given by the identity 

{n{u)\v)n + {ht{ut)\vt)t - ([e(w)l + er(wr)H(r)|/i)r ds = 0, 

which at equilibrium yields 

{ii*\'v)n + {ht*W)t + u^{l*\h)T = 0, (2.6) 

where k, := kt* '■= '«r(w*) and 

— —{i{u,) + ;rK)H(r,)} = iv^'K)! + a'{u,)n{T,). (2.7) 

On the other hand, a straightforward calculation yields with z — (v, v-p, h) 

-{'D^z\z) = — [{i<i*v\v)n + (Kr*i'r|i'r)r - (T*w*('H'(r*)/i|/i)r] , (2.8) 

where kt* = uriu*) and ct* = a{u*). As k* and kt* are positive, we see that the 
form (Vzlz) is negative semi-definite as soon as 'H'(r*) is negative semi-definite. 
We have 

where A* denotes the Laplace-Beltrami operator on and i?* means the radius 
of an equilibrium sphere. To derive necessary conditions for an equilibrium e* to 
be a local maximum of entropy, we consider two cases. 

1. Suppose that T^, is not connected, i.e. is a finite union of spheres F^'. Set 
V = vr — 0, and let h — h^he constant on Fj with — 0. Then the constraint 

(|2.6p holds, and with w„ the area of the unit sphere in R" 

m 

{Vz\z) - (a,^/,)((n- ^/i^ > 0, 

fe=l 
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hence T) cannot be negative semi-dcfinitc in this case, as cr* > by (|2.3p . Thus if 
e* is an equilibrium with maximal total entropy, then must be connected, and 
hence both phases are connected. 

2. Assume that F* is connected. With h = — (K*|l)n — Kr*|r*|, v = vv = m*^*|F*| 
we see that V negative semi-definite on the kernel of E'(e*) implies the condition 

■= ^^iiw:\ ^ ' ^'-'^ 

We will see below that connectedness of F* and the strong stability condition < 1 
are sufficient for stability of the equilibrium e* . 

We point out that the quantity defined in (|2.9I) coincides with the analog 
quantity in [181 Defintion (1.11)] in case Kr* = and a — constant. (Note that 

= l{uj,)/uj, in this case, which differs from the definition of in |18)). 

(f) Summarizing, we have shown 

• The total energy is constant along smooth solutions of (jl.6p . 

• The negative total entropy is a strict Ljapunov functional for (|1.6p . 

• The equilibria of (|1.6I) are precisely the critical points of the entropy func- 
tional with prescribed energy. 

• If the entropy functional with prescribed energy has a local maximum at 

= (M,t,ur*jr*) then F^ is connected. 

• If F* is connected, a necessary condition for a critical point (u*,ur*,r*) 
to be a local maximum of the entropy functional with prescribed energy 
is inequality (I2.9p . 

(g) Now let us look at the energy of an equilibrium as a function of temperature. 
Suppose we have an equilibrium (it, F) at a given energy level Eq, and assume that 
F consists of m disjoint spheres of radius R contained in VL. Then 

< i? < Rm '■— sup{i? > : f2 contains m disjoint ball of radius i?}, 

and with Lp{u) ;= ["(/"(u)] we have 

= ipiu) + cr(u)-H(F) = ^{u) - (n - \)a{u)/R, 

and hence R = R{u) = [n — \)a{u) / ipiu). Further we have 



Ee(u) := E(7i,F)= j e{u) dx + J er{u)ds 



^e2iu)\n\~\k\Hu)]+eliu)\T\ 
= e2iu)\n\ - (mcj„/n)i?(u)"|e(u)l +mcj„i?(u)"-ier(u) 



£2{u)\il\ + Cn,r. 



cr(M)" d a{uY 

u- 



where c„_m = muJn{n — l)"^^/n. Thus we obtain for the total energy of an equi- 
librium 

Ee{u)=5{u)-uS'{u), S{u) = + cn,^ /^^ (2.10) 
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Consequently, the equilibrium temperature for an equilibrium, where F consists of 
m components, is the solution of the scalar problem 

Eq = Eg{u) = 5{u) — u5'{u), < u < Uc, < a{u)/ip{u) < Rm/{n — 1). 

Let us look at the derivative of the function Ee(u). A simple calculation yields 

E.'^{u) = -uS"{u) = -\n\u^'^{u) - Cu,mU^ [n {a/ipr-' a'-{n- l)(a/^)>'] 
= \n\n2{u) - Cn,mu[n{a/^r'^a" - (n - 1)(^/<^)'V"] 

- c„,™n(n - l){a/^r-^u[{a'f/a - 2aV7^ + ^i^'f/v^] 
= \D.\k2{u) + |r|ACr(w) - \^i\{<u)\ - {R^\T\/{Ti-l)a)u[^' - a'^/af 
= [(«(u)|l)o + |r|«rW]-(i?'H|r|/(n-l)a(u))7.[I^'(«)l+a'(u)7^(r)]l 

Therefore the stability condition C,{u) < 1 is equivalent to Eg(u) < 0, an alternative 
interpretation to the one obtained above. 

3. Transformation to a Fixed Interface 

Let c M" be a bounded domain with boundary dO. of class C^, and suppose 
r C is a closed hypersurface of class C^, i.e. a C^-manifold which is the boundary 
of a bounded domain ili C Q.. We then set 0,2 = \ Jli. Note that while ^2 
is connected, fii may be disconnected. However, Qi consists of finitely many 
components only, as dV,i = F by assumption is a manifold, at least of class C^. 
Recall that the second order bundle of F is given by 

Note that the Weingarten map Lr (also called the shape operator, or the second 
fundamental tensor) is defined by 

Lr{p) = -VtMp), per, 

where Vr denotes the surface gradient on F. The eigenvalues Kj{p) of Lr{p) are 
the principal curvatures of F at p € F, and we have |ir(p)| = maxj |/tj(p)|. The 
curvature Hr ip) is defined by 

n-l 

^r{p) = ^ Kj{p) = trLr(p) = -divri^{p), 

i=i 

where divr means surface divergence. Recall also that the Hausdorff distance dn 
between the two closed subsets A,B c is defined by 

dH{A,B) := max{ sup dist(a, B), supdist(6, A)} . 

Then we may approximate T by a real analytic hypersurface S (or merely S € C^), 
in the sense that the Hausdorff distance of the second order bundles of F and S 
is as small as we want. More precisely, for each 77 > there is a real analytic 



ON STEFAN PROBLEMS WITH VARIABLE SURFACE ENERGY 



11 



closed hypersurface such that (ii/(A/'^I], AA^F) < 77. If 77 > is small enough, then 
S bounds a domain Slf with Of c O, and we set O2 = O \ Of . 

It is well known that such a hypersurface S admits a tubular neighborhood, 
which means that there is a > such that the map 

A : E X i-a,a) M" 

A(p,r) := p + rvsip) 

is a diffeomorphism from S x {—a, a) onto TZ{A). The inverse 

A"^ :7e(A) h-)- S X {-a, a) 

of this map is conveniently decomposed as 

A-\x) = (nE(a;), dj:{x)), x e n{K). 

Here H^{x) means the nonlinear orthogonal projection of x to E and d-s{x) the 
signed distance from x to E; so |(iE(a;)| = dist(a;, E) and ds{x) < iff x S flf. In 
particular we have TZ{A) = {x e M" : dist(a;, E) < a}. 

On the one hand, a is determined by the curvatures of E, i.e. we must have 

< a < min {l/|Kj(p)| : j = 1, . . . ,n — 1, p G E}, 

where Kj{p) mean thc^ principal curvatures of E at p G E. But on the other hand, 
a is also connected to the topology of E, which can be expressed as follows. Since 
E is a compact (smooth) manifold of dimension n — 1 it satisfies a (interior and 
exterior) ball condition, which means that there is a radius rs > such that for 
each point p € E there are xj € O?, j = 1,2, such that Brs{xj) C Oj, and 
Brs{xj) n E = {p}. Choosing rs maximal, we then must also have a < r-^. In the 
sequel we fix 

a = ;^min{rE, ] , j = 1, . . . , n - 1, p e eI . 

2 i \Kj{p)\ J 

For later use we note that the derivatives of Ilj:{x) and dj:{x) are given by 

Vdj:ix) = MUj:{x)), Il'^{x) = Mo{dj:{x),U{x))Ps{Uj:{x)) 

for < a, where -Pe;(p) = I—i^i:{p)®i^Y:{p) denotes the orthogonal projection 

onto the tangent space TpE of E at p G E, and 

Mo{r){p) = {I - rLj:{p))-\ (r,p) € (-a, a) x E. (3.1) 

Note that 

|Mo(r)(p)| < l/(l-r|Ls(p)|) < 2, for aU (r,p) e (-a, a) x E. 

Setting r = r(t), we may use the map A to parameterize the unknown free bound- 
ary r{t) over E by means of a height function h{t,p) via 

r{t) = {p + h{t,p)v^{p) : p e E, t > 0}, 

at least for small \h\^. Extend this diffeomorphism to all of O by means of 

'Eh{t,x) =x + x{dT.{x) / a)h{t,IiT.{x))i>Y.iJlY,{x)) =:x + ^h{t, x). 
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Here x denotes a suitable cut-off function. More precisely, x G T^^)i < x < Ij 
X(r) = 1 for \r\ < 1/3, and xW = for \r\ > 2/3. Note that Eh{t,x) = x for 
|d(a;)| > 2a/3, and 

S'f^^{t,x) — X — h(t,x)i>-s{x), a; e E, 

for |/i|oo sufficiently small. 
Setting 

v{t,x) = u{t,'E.p{t,x)) or u{t,x) — v{t,'E.p^{t,x)) 

we have this way transformed the time varying regions n\r(t) to the fixed domain 
ri \ E. This is the direct mapping method, also called Hanzawa transformation. 
By means of this transformation, we obtain the following transformed problem. 

K{v)dtv + A{v, p)v = k{v)TZ{p)v in \ E 

di,v — on dQ 

|wj = 0, vp — V on E 

^ivr)j+<yivr)nip)-jivr)[3{p)dtp = on E 

KT{vr)dtVr +C{vr,p)vr + B{v,p)v = 

= ~{1{vt) + It{vt)H{p) - i{vT)P{p)dtp}P{p)dtp on E 

w(0) = Wo, /5(0) = po- 

(3.2) 

Here A{v, p), B{v, p) and C{vr, p) denote the transformed versions of the operators 
of — div(dV), — and — divr(drVr), respectively. Moreover, Hip) means the 
mean curvature of F, l3{p) = (i/j^lvrip)), the term f3(p)dtp represents the normal 
velocity V, and 

TZ{p)v = dtv - dtu o Sp. 

The system p. 21) is a quasi-linear parabolic problem on the domain 17 with fixed 
interface E C i7 with dynamic boundary conditions. 

To elaborate on the structure of this problem in more detail, we calculate 

DEp = I + DCp, [DEp]-' =I-[I + D^p]-'D^p / - A'hip)'^ . 

where D deontes the derivative with respect to the space variables. Hence — 
for |(iE(a;)| > 2a/3 and 

DCp{t,x) = -x'(ds(a;)/a)p(i,Hs(x))i.s(Hs(a;)) ® i^s(Hs(x)) 
a 

+ x(ds {x)/a) (Hs (x)) ® Mo{ds (x)) Vsp(t, Hs(a;))] 
- x(ds {x)/a)p{t, Hs (Hs (a:))Mo (ds (a:))PE (He (a:)) 

for < |(iE(a;)| < 2a/3. In particular, for x G E we have 

D^p{t,x) = i^s{x) (g) VEp(i,a;) - p{t,x)L^{x)P^{x), 

and 

[D^pVit, x) = Vsp(t, x) (g) zys(a;) - p(t, x)Lj:{x), 
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since L-s:{x) is symmetric and has range in T^T,. Therefore, [I + DS^p] is boundedly 
invertible, if p and Vsp are sufficiently small, and 

Employing this notation we obtain 

Vu o Ep = {[DE;')]' o Ep)Vv = [DEp]-''^Vv =: (/ - Mi(p))Vu, 
and for a vector field q = q oEp 

(V|g) o Ep = (([i^S; Y o Ep)V\q) = i[DEp]-''-'V\q) = ((/ - M,{p))V\q). 
Further we have 

dtu oEp = dtv - (Vm o Ep\dtEp) = dtv - ((DS^ ^ o Ep}Vv\dtEp) 

= dtV - {[DEp]-'^^Vv\dt^p) = dtV - {Vv\{I - Mj{p))dt^p), 

hence 

n{p)v = {Vv\{i-Mj{p))dtCp). 

The normal time derivative transforms as 

dt,nUr oEp = dtvr + (Vsfrlz^)!^ = dtvr, 

as Vs^r is perpendicular to v-^. 

With the Weingarten tensor is = — V^i^e we obtain 

^{p) = P{p){vi: - oi{p)), a{p) = Mq{p)Vj:P, 

Mo(p) = (/ - pLj:)-\ Pip) = (1 + |a(p)|2)-^/^ 

and 

V = {dtEp\iyr) = {yT,\vT{p))dtp = P{p)dtp. 
For the mean curvature T-L{p) we have 

^(p) = /?(p){tr[Mo(p)(is + Vs«(p))]-/?'(p)(Mo(p)a(p)|[Vsa(p)]a(p))}, 

an expression involving second order derivatives of p only linearly. More precisely, 

H{p) = P{p)Qip)-^lp+m^{p), 

Gip) = Mlip) - p^{p)Moip)V^p (g> Mo(p)Vsp. 

Note that /3 as well as and G only depend on p and Vi;p. The linearization of 
the curvature at p = is given by 

W'(0) =tri| + AE. 

Here As denotes the Laplace-Beltrami operator on S. B{v, p) becomes 

B{v,p)v = -ld{u)d,ul o Ep = -m^W - M^{p))Wv\\vt) 
= -p{p){{d{v){I - M^{p))Vv\\u^ - a{p)) 
= -l3{p)ld{v)d,,vl+l3{p){ld{v)Vvl\{I - M,{p)fa{p)), 
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since MJ{p)vy: = 0, and 

Aiv, p)v = - div(d(u)VM) oEp = -{{I - Mi(p))V|d(w)(/ - Mi(p))V'y) 
= - d{v)Av + d(v)[Mi{p) + Mj{p) - Mi{p)Mj{p)] : W^v 
- d'{v)\il - Mi(p))Vt;p + d{v){{I - Mi(p)) : VMi(p)|Vi;). 

We recall that for matrices A,B £ K"^", A : B ^ T,lj=iaijbij = ti {AB'^) 
denotes their inner product. The pull back of Vr is given by 

Vr<y3 oEp^ Pr{p)MQ{p)V^ip, 

where 

Pr{p) = I - vt{p) ® vt{p). 
This implies for C{vy, p)vr the relation 

Civr,p)vr = -tr{Pr(p)Mo(p)VE(dr(fr)/'r(p)Mo(p)Vswr)}. 

It is easy to see that the leading part of A(v, p)v is —d{v)Av, while that of B{v, p)v 
is — [}{p)\d{v)d^^v\, and the leading part of C(wr, p)vt turns out to be — dr(wr)As- 
This follows from A/o(0) = 1, Pr(0) Pe, A/i(0) = and a(0) = 0; recall that we 
may assume p small in the C^-norm. It is important to recognize the quasilinear 
structure of 



(4.1) 



4. Linearization at Equilibria 

The full linearization at an equilibrium (u*,Mr*,r*) with ur* — T* = UfeEfe 
a finite union of disjoint spheres contained in Q, and with radius i?* > given by 
R^ = {n— 1)(t(u*)/|V'(m*)], reads 

K^,dtv — d^:Av — K^,f in \ 
d,yV — on dil 

{vj ^ 0,vr — V on T^, 

KT*dtVT — dr* A*wr — |d*c?i/ti] + hu*dtp ~ Kr*/r on F* 
hvr — <J*A^,p — j^,dtp — g on F* 

v{0) = vq, p{0) = po. 

Here 

K* — > 0, Kr* — Kr(u*) > 0, = d{u^,) > 0, 

dr* ~ rfr(u*) > 0, tr* = cr(M*) > 0, 7* = l{u*) > 0, 

and as in 

= |V''(u*)l + o-'(u*)H(F*) = — (t'(u*)(p(u*)/o-(w*) = (t(u*)A'(m*), 

and 

^* = -(-^ + ^*)' 
where A, denotes the Laplace- Beltrami operator on F*. 



ON STEFAN PROBLEMS WITH VARIABLE SURFACE ENERGY 



15 



4.1. Maximal Regularity. We begin with the case 7* > 0, which is the simpler 
one. Define the operator L in 

with 

Xi w^{n \ r,) X VK2+'^(r,) x w^+'{t,) 

by means of 

D{L)^{{v,vr,p)CzXi: H = 0, ur = w on T*, d^v^QoudVt), 

(-d,/K,)A 

L= — |(d*/Kr*)<9i/] (/^u*/7* — dr*A,)/Kr* — (/»w*cr*/7*Kr*)^* 
-(/,/7*) (cr*/7,)A* 

In case 7* > 0, problem (14.11) is equivalent to the Cauchy problem 

z + Lz = (/, /r - {Uu^/"f*KT*)g,g), z(Q) = zq, 

where z = (D,wr,p) and zq = (wq, uoIfoj Po)- The main result on problem (j4.ip for 
7* > is the following. 

Theorem 4.1. Let 1 < p < 00, 7* > 0, and 

-1/p < r < 1 - 1/p, r<s<r + 2. 
Then for each finite interval J — [0, a], there is a unique solution 

{v,vr,p)eE{J) -.^ HliJ:Xo)nLp{J;Xi) 
of (j4.ip if and only if the data {f,fTig) and {vo,vro,po) satisfy 
(/,/r,g)eF(J) = Lp(J;Xo), 

{vo,vro,Po) e W^-'/Pin\T^,) X W^+^-'^p{T^,) x W^+^~'^p(T,) 

and the compatibility conditions 

[uol = 0, UFO = '^0 on r*, di,v = on dil. 

The operator —L defined above generates an analytic Co-semigroup in Xq with 
maximal regularity of type Lp . 

Proof. Looking at the entries of L we see that L : Xi Xq is bounded provided 
r < 1 — 1/p, r < s, and s < r + 2. The compatibility condition v-p = v^^^ implies 
r + 2 > 2— l/p. This explains the constraints on the parameters r and s. To obtain 
maximal Lp-regularity, we first consider the case s > r. Then L is lower triangular 
up to perturbation. So we may solve the problem for {v,v-r) with maximal Lp- 
regularity (cf. [1] for the one-phase case) first and then that for p. In the other 
case we have r = s. Then the second term in the third line in the definition of L is 
of lower order, hence p decouples from (v,v-p). This way we also obtain maximal 
Lp-regularity. Since the Cauchy problem for L has maximal Lp-regularity, we can 
now infer from (14t Proposition 1.2] that — L generates an analytic Co-semigroup 
in Xq. □ 
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We note that if = and 7, = then the hnear problem (|4.1|) is not well-posed. 
In fact, in this case the linear Gibbs- Thomson relation reads 

which is not well-posed as the kernel of A* is non-trivial and A^, is not surjective. 

Now we consider the case Z* ^ and 7* — 0. For the solution space we fix again 
r, s e M with r < s < r + 2, —1/p < r < 1 — 1/p, and consider 

{v,vr,p) eEiJ) = HliJ,Xo)nLpiJ;Xi). 

Then by trace theory the space of data becomes 

ifjr,g) e Fo(J) Lp(J;Lp(r!)) x Lp{J;W;{T.,)) 

X [H^{J; W;-'{T,) n Lp{J; W;{T,))], 

and the space of initial values will be 

K,«ro,Po) e w^-'/p{n\T,) X iy;+2-2/P(r,) x w;+'-^^piT,) 

with compatibilities 

|t)o] = 0, vro = vq, Z*vro - a'*^*Po = 5(0) on T*, d^v = on dQ. 

To obtain maximal Lp-regularity, we replace vr by the Gibbs-Thomson relation, 
which for 7* = is an elliptic equation. We obtain vr = (cr*//*)A*p -t- g/l*. 
Inserting this expression into the energy balance on the surface yields 

-I- {Kr*(7^/h)A^)dtp ~ dr*A^vr - Id^d^vj = Kr*(/r - dtg/h)- (4.2) 
Moreover, we obtain 

— {hu^dr*/Kr*)A^p + {dY*/l*)A^g. 

Now we assume that 

u^;2^2 _^ (n- l)Kr*cr*, (4.3) 
which is equivalent to invertibility of the operator Aq :— l^u^ + {KY*cr*/h)A^. 
Applying its inverse to ()4.2p we arrive at the following equation for p: 

dtp - (dr*/Kr*)A*p + ^{(u*/*dr*/Kr*)A*/3 - = g, (4.4) 

with 

g = Aq ^{Kr*/r - {{nr*/l*)dtg - {dr*/h)A^g)}. 
Solving equation (|4.2p for dtp we obtain for vr'- 

Kr*dtvr ~ dr*A*wr — + Z*u*Aq ^{dr*A*i;r + = /r- (4.5) 

where 

/r = Kr*{/r - ^*u*Ao - dtg/h)}. 
Then by the regularity of /r and g and with r < s < r + 2 we see that 
/r e LpiJ; W;iT,)), g e Lp(J; M/;(F,)). 
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So the linear problem ()4.1|) can be recast as an evolution equation in Xq as 

Z + LqZ 

with Lq = Lqq + Lqi defined by 
D{Loj) = {{v,vr,p) e Xi : 



= {fjr,g), z(0) = zo, 
|d] = 0, fr = f on r*, d^v = on dfl}, 



and 



^00 



and 



(-d*/K*)A 
-I(d*/«:r*)9.] 







-((ir*/Kr*)A+ 











-((ir*/Kr*)A, 





{hu^c/ Kr*)AQ^ld^d,y] (hu^dr*/ Kr*)A^^ 

(u*Z*dr*/Kr*)^(7^A* 

Looking at Lq wc first note that Lqi is a lower order perturbation of Lqq. The latter 
is lower triangular, and the problem for {v, vr) as above has maximal Lp-regularity 
in Xq. As the diagonal entry in the equation for p has maximal Lp-regularity as 
well we may conclude that —Lq generates an analytic Co-semigroup with maximal 
regularity in Xq More precisely, we have the following result. 

Theorem 4.2. Let 1 < p < oo, ^ 0, -l/p < r < I - 1/p, r < s < r + 2, 
h 0, and assume u.J^R'^ ^ Kr*cr*(n — 1). 

Then for each interval J = [0,a], there is a unique solution (v,vr,p) G ^{J) of 
J-^.Jp if and only if the data (/, /r,g) and (t'cfpcPo) satisfy 

(/,/r,.9)eFo(J), 

K,«ro,Po) G w^-'/p{n\T,), xw;+2-2/p(r,) x w;+'-^^p{t.,) 

and the compatibility conditions 

fuoj — 0, vro — I'o, l^,vo — a^,A^,po — g{0) on F*, d^v = on dH.. 

The operator ~Lq defined above generates an analytic Co-semigroup in Xq with 
maximal regularity of type Lp . 

Note that the compatibility condition l^,vo — a^,A^po ~ g{0) allows to recover 
the Gibbs-Thomson relation from the dynamic equations. Indeed, it follows from 
()4.4|) - ()4.5|) that the function w := vr — {{cr*/h)A^p + g/h) satisfies 

KT*dtw — dr*A^w — hu^[dtp - ^{fir*A*ur + + Kr*(/r — dtg/l*)}) 

— —Uu^AQ^dr*A^.w, 

and this shows that w is a solution of the parabolic equation 

Kr*dtw ~ dr^A^vu + l^u^dr^A^^A^w = 0, w{0) = on F*. (4.6) 

As w = is the unique solution of ()4.6p we conclude that the Gibbs-Thomson 
relation is satisfied. 
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4.2. The Eigenvalue Problem. By compact embedding, the spectrmii of L con- 
sists only of countably many discrete eigenvalues of finite multiplicity and is inde- 
pendent of p. Therefore it is enough to consider the case p = 2. In the following, 
we will use the notation 



{9\h)T, ■■= {g\h)L2{T.)-- 



IV dx, It, w e L2(ri), 
gh ds, g,h e i2(r*), 



r. 



/ I ^l/2 



for the L2 inner product in Q and F*, respectively. Moreover, we set \v\n 
1/2 

and l^lr. ~ (sl.?)/, ■ The eigenvalue problem reads as follows: 

K*Aw — d*Au = in \ 
d,yV = on dfl 
Ivj = on F, (4.7) 

hv — a^,A^,p — 'j^.Xp ~ on F* 
Kr*Xv — dr*A^,v — ld^,d^vj + hu^,Xp — on F*. 

Let A 7^ be an eigenvalue with eigenfunction {v,p) 7^ 0. Then (|4.7p yields 
= \\y/^v\l - {d*Av\v)n = X\y/^v\l + \y^Wv\l + {ld,d^vj\v)r,. 
On the other hand, we have on the interface 

= Kr*A|w|r. — dr*(Arw|w)r. — (N*c^!^«ll^')r. + XuJr{p\v)r, 
==AKr*|w|r. +c^r*|Vrf|r. ^ {ld*d^v]\v)r, -I- Au*/*(p|w)r. . 
Adding these identities we obtain 

= X\^/K^v\fi + \^/d^Vv\fi + AKr*|i'|r. + '^r*|Vrw|r. + Au,^*(p|u)r. , 
hence employing the Gibbs- Thomson law this results into the relation 

+ Xutar{A^,p\p)r, -I- 7*u*|A|^|p|r, =0. 

Since A* is selfadjoint in L2(r*), this identity shows that all eigenvalues of L are 
real. Decomposing v = vq + v, vr = ur.o + vr, p = Pa + P, with the normalizations 
(K*|wo)n — (wr,o|l)r. ~ (po|l)r. = 0, this identity can be rewritten as 

X{\y/K^vo\n + Kr*|wr,olr. + (^*u*{A.^.po\po)r, + Au*7*|polr. } 



+ A 



(k,|1)u^ + Kr*|r*|wp — o-,u» 



1 



F,|p^ + 7*M*A|r*|/9^ 



= 0. 



In case F* is connected, A^ is positive semi-definite on functions with mean zero, 
and hence the bracket determines whether there are positive eigenvalues. Taking 
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the mean in (|4.7p we obtain 

{K*\l)nv + Kr*\^*\vr + /*M*|r*|p — 0. 
Hence minimizing the function 

with respect to the constraint we see that there are no positive eigenvalues provided 
the stabihty condition < 1 is satisfied. 

If r, = Ui</<m consists of m > 1 spheres F'^ of equal radius, then 

N{L) = span|( ^*^"^^^^"% -/0, (0, YI), . . . , (0, F^) : 1 < Z < m| , (4.8) 

where the functions denote the spherical harmonics of degree one on F'^ (and 
Y^' = on {Ji-^i^D, normahzed by (YJIYDyi^ = Sjk- N{L) is isomorphic to the 
tangent space of £ at (uh.,F*) S £, as was shown in [18l Theorem 4.5.(vii)]. 
We can now state the main result on linear stability. 

Theorem 4.3. Let cr, > 0, 7, >0,h^ 0, 

ry* := {n — l)(T*Kr*/u*Z^i?^ ^1 in case — 0, 
and assume that the interface F, consists 0/ m > 1 components. Let 

_ [n - 1)o-»[(k^|1)o + Kr*|r^|] 
^* " uJimVA ' 

and let the equilibrium energy Eg he defined as in (j2.10p . Then 

(i) E'M = (C* - l)uJlRl\T,\/{n~-l)a^. 

(ii) Q is a an eigenvalue of L with geometric multiplicity (mn + 1). 

(iii) is semi-simple if C* 7^ 1 • 

(iv) //F* is connected and ^* < 1, or if 77, > 1 and 7* = 0, then all eigenvalues 
of —L are negative, except for the eigenvalue 0. 

(v) // > 1, and 77* < 1 m case 7* = 0, then there are precisely m positive 
eigenvalues of —L, where m denotes the number of equilibrium spheres. 

(vi) //C* ^ 1, ctnd 77* < 1 m case 7* = then —L has precisely m — 1 positive 
eigenvalues. 

(vii) N(L) is isomorphic to the tangent space T(^u,,r»}^ '^f ^ '^^ (m*,F*) G £. 

Remarks 4.4. (a) Formally, the result is also true if = and 7* > 0. In 
this case Eg(u*) = («;*|l)f2 + Kr*|r*| > and = 00, hence the equilibrium is 
unstable. If in addition 7* — 0, then the problem is not well-posed. 

(b) Note that C* does neither depend on the diffusivities d*, dr*, nor on the coef- 
ficient of undercooling 7*. 

(c) It is shown in that in case C* = 1 and F* connected, the eigenvalue is no 
longer semi-simple: its algebraic multiplicity rises by 1 to (n + 2). 
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(d) It is remarkable that in case kinetic undercooling is absent, large surface heat 
capacity, i.e. 77* > 1, stabilizes the system, even in such a way that multiple spheres 
are stable, in contrast to the case ry* < 1. 

(e) We can show that, in case 7* = 0, if 77* increases to 1 then all positive eigen- 
values go to CX). 

We recall a result on the Dirichlet-to-Neumann operator D\, A > which is 
defined as follows. Let g G i/2^^(r*) be given. Solve the elliptic transmission 
problem 

' K^Aw — d*Aw = in \ r*, 
duW = on 9fi, 

M=0 onr„ ^^'^^ 

w ^ g on r*, 

and define Dxg = e H^^^{r^). 

Lemma 4.5. The Dirichlet-to-Neumann operator D\ has the following well-known 
properties. 

(a) {Dxg\g)T, = X\^l'^w\l + \dl'^Vw\l, for all g e hI'\t,); 

(b) \Dxg\r, < C[X^/^\g\r. + Iglniir,)], for all g e i?2'^'(r,) and A > 1; 

(c) (i:'A5|5)r. > cAi/2|g|2^, for all g e h!^^^{T,) and A > 1. 

In particular, D\ extends to a self adjoint positive definite linear operator in Li (F*) 
with domain H\{Vtf). 

4.3. Proof of Theorem 14.31 For the case that Kr* — dr* = this result is 
proved in [T5^. Assertion (i) follows from the considerations in part (g) of the 
introduction. Assertions (i), (iii), and (vii) only involve the kernel of L and the 
manifold of equilibria. Since both are the same as in the case Kr* = c?r* = 0, the 
proofs of (i), (iii) and (vii) given in |18j remain valid in the more general situation 
considered here. The first part of assertion (iv) has been proved above, it thus 
remains to prove the assertions in (v) and (vi), and the second part of (iv). 

If the stability condition C* < 1 does not hold or if F* is disconnected, then there 
is always a positive eigenvalue. It is a delicate task to prove this. The principal 
idea to attack this problem is as follows: suppose A > is an eigenvalue, and that 
p is known; solve the resolvent diffusion problem 

K*Aw — Aw = in \ F* 

dyV = on 9ri 

\v\ - on F, ^^'^^^ 
V = vy on F:,, 

to get — =: D\vr- Next we solve the resolvent surface diffusion problem 
XKr*vr — c?r*A*wr + D\vr — h, 



ON STEFAN PROBLEMS WITH VARIABLE SURFACE ENERGY 



21 



to the result 

vr = Txh := (A^r* — rfr*A* + Dx)^^h. 
Setting h = —XuJ^,p this implies with the linearized Gibbs- Thomson law the 
equation 

[(^*u*)ArA+7*A]p + cr*Ap = 0. (4.11) 
A > is an eigenvalue of — L if and only if (14. lip admits a nontrivial solution. We 
consider this problem in L2(r*). Then A.^ is selfadjoint in L2(r») and 



for each p E D{A^,) — i?|(r,). Moreover, since A^, has compact resolvent, the 
operator 

Bx := [{llu.)XTx + 7*A] + a,A, (4.12) 
has compact resolvent as well, for each A > 0. Therefore the spectrum of Bx 
consists only of eigenvalues which, in addition, are real. We intend to prove that 
in case either ^^, is disconnected or the stability condition does not hold, Bxg has 
as an eigenvalue, for some Aq > 0. This has been achieved in [18^ in the simpler 
case where Kr* = dp* = 0, in which case Tx is the Neumann-to-Dirichlet operator 
for ()4.10|) . Here we try to use similar ideas as in flF, namely we investigate Bx 
for small and for large values of A. However, in the situation of this paper this will 
be more involved. 

For this purpose we need more information about Tx- So we first consider the 
problem 

Kf,Xv — d■^,Av = in \ 
— on dft 
Ivj =0 on F, ^^-"^^^ 
AKr*w — dr*AtV — |d*9,yw] — g on F*. 
As we have seen above this problem has a unique solution for each A > 0, denoted 
by V = SxQ- Obviously for A = this problem has a one-dimensional eigenspace 
spanned by the constant function e = 1. The problem is solvable if and only if 
the mean value of g is zero, i.e. if g £ i2,o(r*). This implies by compactness that 
Sxg — Sog as well as Tx — Tgg as A 0+, whenever g has mean zero, where S^g 
means the unique solution of (|4.13p for A = with mean zero. 

(a) Suppose that F* is disconnected. If the interface F* consists of m components 
FJ, k = 1, m, we set e/c = 1 on F^ and zero elsewhere. Let p = 0^6^ ^ with 
X]fc ^fe = 0> hence Qop — p, where Qo is the canonical projection onto L2.o(F,) in 
i2(r*), QoP ■=P- (p|e)r./|F,|. Then 

lim XTxp — lim XTxQoP = 0, 
since TxQq is bounded as A — > 0. This implies 

\im {Bxp\p)r, = - {{n - l)a,/Rl) ^ \Vl\al < 0. 

k 
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(4.14) 



Therefore B\ is not positive semi-definite for small A. 
(b) Suppose next that is connected. Consider p = e. Then we have 
iBxe\e)r,^uJlXiTxe\e)r,+XjMl,-iin~l)a,/Rl)\e\l^. 
We compute the limit limA^o A(7Ae|e)r, as follows. First solve the problem 

— d*Au — K*ao in \ 
d„v — on dfl 

Iv] =0 on F* 

-dr*Arw — |rf*9i/w] = e — Kr*ao on F*, 

where ao = |r*|/[(K*|l)o -l-Kr*|r*|], which is solvable since the necessary compati- 
bility condition holds. Let vq denote the solution which satisfies the normalization 
condition (re*|t'o)o + Kr*(wo|l)r, = 0. Then v\ := S\e — vq — cq/A satisfies the 
problem 

K^,Xvx — d^,Avx = —k*Xvq in \ F* 
d,^v = on dfl 

Ivx] = on F, 

^ KT*\v\ — (ir*A*i;A — |d*9iyfA] — — AKr*fo on F*. 

By the normalization (K*|wo)f2 + '*r*(i'o|l)r, = we see that the compatibility 
condition for (|4.13p holds for each A > 0, and so we conclude that v\ is bounded 
in W|(ri \ F*) as A — ?► 0, it even converges to 0. Hence we have 

lim XTxe = lim [{Xv\ + Auo)|r. + oo] = ao. 



(4.15) 



This then implies 

lim(BAe|e)r, = l^u^ 



{n — l)o-, |F* 



<o, 



(«:*|l)a + Kr*|F*| 

if the stability condition does not hold, i.e. if C* > 1- Therefore also in this case 
B\ is not positive semi-definite for small A > 0. 

(c) Next we consider the behavior of {Bxp\p)r, as A — > oo. We intend to show 
that Bx is positive definite for large A. We have 

ATa = A(Kr*A - dr*Ar + I^a)"^ ^ 1/Kr* for A ^ oo, 

as Dx is of lower order, by (b) of Lemma 4.5. This implies for a given g E D(A^,) 

{B\g\g)r. = llu.A{Txg\g)r , + (T^{A.^.g\g)T* +7*A|g|r. 



> (7*A 
^ (7*A- 



(n — l)cr* 

m 

{n - l)cr, 

m 



)\g\l,+i:uA{Txg\g)r, 

2„ 
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as A — > oo. We have thus shown that B\ is positive definite if 7* > and 
A > (n — l)(T*/7*i?^, or if 

7* = and IIu^/kt* > (n - l)(T*/i?^ (4.16) 

In particular, for 7* = and small the latter condition condition will be violated, 
in general. 

(d) In summary, concentrating on the cases 7* > or (|4.16p , we have shown that 
B\ is not positive semi-definite for small A > if either r* is not connected or 
the stability condition does not hold, and B\ is always positive definite for large 
A. Let 

Ao = sup{A > : Bfj^ is not positive semi-definite for each /i G (0, A]}. 

Since Bx has compact resolvent, Bx has a negative eigenvalue for each A < Aq. 
This implies that is an eigenvalue of Bxo, thereby proving that ~L admits the 
positive eigenvalue Aq. 

Moreover, we have also shown that 

Bop = lim [llu^XTxp + 7*Ap + a^A^p] = - — * * * Pop + cr^A^p, 

where P^p {I ~ Qo)p — (p|e)r,/|r*|. Therefore, Bq has the eigenvalue 
uJl\T*\ (n - l)a^ _ uJl\T^\ ^ ^ ^ 

with eigenfunction e, and in case m > 1 it also has the eigenvalue —{n — l)(T*/i?J 
with precisely to — 1 linearly independent eigenfunctions of the form ak^k with 
Efcflfe = 0- 

As A varies from to Ao, all the negative eigenvalues of Bq identified above 
will eventually have to cross along the real axis. At each of these occasions, —L 
will inherit at least one positive eigenvalue, which will then remain positive. This 
implies that —L has exactly to positive eigenvalues if the stability condition does 
not hold, and to — 1 otherwise. This covers the case 7, > as well as (|4.16p . 

(e) To cover the remaining we assume 7, = and Kr*{n~ 1)/Rl > w*?^/(7, =: (5*. 
Suppose A > is an eigenvalue of io- Then there is p ^ such that 

(Akf* — dr* A» + D\)A^,p + \6^,p = 0. 

Multiplying this equation in L2(r*) by A^p and integrating by parts one obtains 
the identity 

AKr,|A,/9|r^ -fdr*|Vr.A,p|r^ + {DxA^p\A^p)r, + X6*{A^p\p)r, =0. 
As Dx is positive definite in ^2(1"*) this equation implies 
AKr*|A*p|p^ + X6t,{At,p\p)r, < 0. 
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Let P denote the projection onto the kernel A/'(A*) and Q — I — P. Since P, Q 
commute with this imphes 

AKr*|A*Qp|r. + XKr*\A^Pp\l^ + XS^{A^Pp\Pp)r, < 0, 

as A^ is positive scmi-dcfinitc on TZ{Q) = TZ{A^). Now A^P = — ((n — 1)/R^)P 
and 



> \nTM*Pp\T, + \5M*Pp\Pp)T. = A-^ ^r*^^ - 6* \Pp\r, > 0, 



n — 1 



hence Pp — and A^Qp = 0. This imphes A^.p — and therefore p = as 
5» > 0. This shows that there are no positive eigenvalues of Lq in case 7* = and 
Kr*{n — 1)/Rl > uJl/a^. This completes the proof. 



5. The Semiflow in Presence of Kinetic Undercooling 

In this section we assume throughout 7(5) > for all < s < Uc, i.e. kinetic 
undercooling is present at the relevant temperature range. In this case we may 
apply the results in [T7] and [TU], resulting in a rather complete analysis of the 
problem. 



5.1. Local Well-Posedness. To prove local well-posedness we employ the direct 
mapping method as introduced in Section 3. As base space we use 

and we set 

X, = {{v,vr,p) e H^p{n\j:) x w^-'/p{e) x w^-'^p{^) : 

H =0, vr = vi^, =0}. 
The trace space X^ then becomes for p > n + 2 

X^ = {iv,vr,p) G W^~^/Pin\J:) X M^p2-3/p(5^) X ^^^^-^/^(E) : 

Ivj = 0, = vi^, d^v\g^^ = 0}, 
and that with the time weight t^~^^, 1 > p > 1/p, 

X^,^ = {iv,vr,p) € W^^-^/Pin\J:) X I^2/.-3/p(j.) X w/2;.+i-3/p(5.) . 

|w] =0, vr = v\^, = 0}, 

Note that 

X^^f, ^ BUC^n \ E) X Ci(E) X C2(E), (5.1) 

provided 2p > 1 + {n + 2)/p, which is feasible as p > n + 2. In the sequel, we 
only consider this range of p. We want to rewrite system p.2p abstractly as the 
quasilinear problem in Xq 

i + A{z)z = F{z), z(0) = zo, (5.2) 
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where z = (ti,wr,p) and zq = (wo, '^ro, Po)- Here the quasihnear part A{z) is the 
diagonal matrix operator defined by 



-A{z) — diag 



(dHM«))(A-M2(p):V2) 

2 

S 



(^(«r)/7(«r))a(p) : 



with M2{p) = Afi(p) + Mj^p) - Afi(p)M7(p)- The seniihnear part F[z) is given 

by 

Tl^P> + - A/i(p))Vi;f - d[y){{l - Mi(p)) : VAfi(p)|V«)} ' 

^{SK, p)« - {[l{vv) + ;rK)H(p) - 7(«r)/?(p)5tp]/?(p)ftp + mg} 

'/'(^^r)//3(p)7(t'r) + a(vY)T{p)h{vv) 

where (p(s) = ^'iid 

TO3 = -dr(wr)(/r(/?)Mo(/'))^ V|t;r - /o)wr- 

We note that 7713 depends on wr, Vswr, and on p, Vsp, V|p, but not on V|tir, 
hence is of lower order. Apparently, the first two components of F(z) contain the 
time derivative dtp\ we may replace it by 

atp = {^(«r)+H(p)}//3(p)7(«r), 

to see that it is of lower order as well. 

Now fix a ball B := Bx^ i^{zo, R)) C Xj^^^, where |po|ci(S) < V for some suffi- 
ciently small 77 > 0. Then it is not difficult to verify that 

provided c?i , T/^i , dr , cr, 7 e C'^(0, 00) and dj,Kj,a,j > on (0,Uc), j — l,2,r, and 
provided 2 > 2/i > 1 + rt + 2/p as before. Moreover, as A{z) is diagonal, well- 
known results about elliptic differential operators show that A{z) has the property 
of maximal regularity of type Lp, and also of type Lp^^, for each z € B. In fact, 
for small j] > 0, A(z) is small perturbation of 

A#{z) = diag[ - {div)/K{v))A, -(dr(t'r)/Kr(fr))As, -((T(t;r)/7(t'r))As] ■ 

Therefore we may apply {IQi Theorem 2.1] to obtain local well-posedness of (|5.2p . 
i.e. a unique local solution 

z e HlJiO, a);Xo) n Lp,^((0, a);X,) ^ C([0, a]; X^,^) n a((0, a];X^) 

which depends continuously on the initial value zq G B. The resulting solution 
map [zq I— !• z{t)] defines a local semiflow in X^,f^. 
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5.2. Nonlinear Stability of Equilibria. Let e* — (u*,Mr*,r*) denote an equi- 
librium as in Section 4. In this case we choose S = F* as a reference manifold, and 
as shown in the previous subsection we obtain the abstract quasilinear parabolic 
problem 

z + A{z)z ^ F{z), z{0) = zo, (5.3) 

with Xq, Xi, Xj as above. We set = (u*, ur*, 0). Assuming that 7^ in the 
stability condition, we have shown in Section 4 that the equilibrium is normally 
hyperbolic. Therefore we may apply [171 Theorems 2.1 and 6.1 ] to obtain the 
following result. 

Theorem 5.1. Let p > n + 2. Suppose 7 > on (0,Uc) and the assumptions of 
(|1.5p hold true. As above £ denotes the set of equilibria of (|5.3p . and we fix some 
z^ £ E. Then we have 

(a) // F:, is connected and C,.^, < 1 then z^, is stable in X-y, and there exists S > 
such that the unique solution z(t) of I15.3\) with initial value zq € X^ satisfying 
\zo — z^\j < S exists on M+ and converges at an exponential rate in X^ to some 
Zoo £ £ as t CO. 

(b) IfT^ is disconnected or if (^^ > 1 then is unstable in X^ and even in Xq. 
For each sufficiently small p > there is 6 €z (0, p] such that the solution z(t) of 
(|5.3p with initial value Zq G X^ subject to \zo — 2*|^ < S either satisfies 

(i) distx.^ (^(^o); ^) > P for some finite time to > 0; or 

(ii) z{t) exists on M-j, and converges at exponential rate in X-y to some z^o G £■ 

Remark 5.2. The only equilibria which are excluded from our analysis are those 
with = 1, which means Eg(w*) — 0. These are critical points of the function 
Ee(u) at which a bifurcation may occur. In fact, if such u^, is a maximum or a 
minimum of Eg then two branches of £ meet at u*, a stable and and an unstable 
one, which means that {u■^,, T^) is a turning point in £. 

5.3. The Local Semiflow on the State Manifold. Here we follow the approach 
introduced in |ll] for the two-phase Navier-Stokes problem and in [18] for the two- 
phase Stefan problem, see also [10| for the MuUins-Sekerka problem. 

We denote by M'H?{^) the closed C^-hypersurfaces contained in Vt. It can be 
shown that ^A'H^{^) is a C^-manifold: the charts are the parameterizations over 
a given hypersurface E according to Section 3, and the tangent space consists of 
the normal vector fields on S. We define a metric on M'H^{n) by means of 

d^„2(Ei,I]2) :=rfff(AA2Si,AA2S2), 

where dn denotes the Hausdorff metric on the compact subsets of R" introduced 
in Section 2. This way MT-l?{^) becomes a Banach manifold of class . 

Let dY,{x) denote the signed distance for E as in Section 2. We may then define 

the canonical level function ip^: by means of 

^s{x)^(l){ds{x)), .xeR", 



ON STEFAN PROBLEMS WITH VARIABLE SURFACE ENERGY 



27 



where 

0(s) = sx{s/a) + (1 - x(s/a))sgns, s G M. 

Then it is easy to see that S = (^^^"'^(O), and V(y5s(a;) — z^s(x), for a; G S. Moreover, 
is an eigenvalue of V'^(fY.{x), and the remaining eigenvalues of V'^tf^{x) are the 
principal curvatures of S at a; G S. 

If we consider the subset MH^ {Q,^r) of ^A'h?{yl) which consists of all closed 
hypersurfaces F G MH^iVt) such that F C satisfies a (interior and exterior) ball 
condition with fixed radius r > 0, then the map 

T : Mn^{n, r) C\n), T(F) := ipr, (5.4) 

is an isomorphism of the metric space A^'H^(r2,r) onto T(A^H^(r2, r)) C C^{d). 

Let s — (n — l)/p > 2. Then we define 

W;{n,r) := {F G MH\n,r) : € W;{n)}. (5.5) 

In this case the local charts for F can be chosen of class as well. A subset 
A C Wp{fl,r) is said to be (relatively) compact, if T{A) C Wp{fl) is (relatively) 
compact. 

As an ambient space for the state manifold of (|1.6p we consider the product 
space C{G) x MV?, due to continuity of temperature and curvature. 
We define the state manifold SM for (I1.6P as follows: 

< u < Uc in ri, dijU = on Silj. 

Charts for this manifold are obtained by the charts induced by M'H^{Vl) followed 
by a Hanzawa transformation as in Section 3. Note that there is no need to 
incorporate the dummy variable ur into the definition of the state manifold, as 
Mr = u|r whenever ur appears. 

Applying the result in subsection 5.1 and re-parameterizing the interface re- 
peatedly, we see that (|1.6|) yields a local semiflow on SM. 

Theorem 5.3. Let p > n + 2. Suppose 7 > on (0,Mc) o,nd the assumptions of 
(jl.Sp hold true. 

Then problem m.6]} generates a local semiflow on the state manifold SA4 . Each 
solution (u,T) exists on a maximal time interval [0,i*), where = t*(Mo,ro). 

5.4. Global Existence and Convergence. There are several obstructions to 
global existence for the Stefan problem with variable surface tension (jl.6p : 

• regularity: the norms of u{t) or F(i) become unbounded; 

• well-posedness: the temperature may reach or 

• geometry: the topology of the interface changes; 
or the interface touches the boundary of fJ; 

or the interface contracts to a point. 
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Let (u, r) be a solution in the state manifold SM. By a uniform ball condition we 
mean the existence of a radius tq > such that for each t, at each point x G r(t) 
there exist centers Xi e fli{t) such that Br^ixi) C Vti and T{t) n Brg{xi) — {x}, 
i = 1,2. Note that this condition bounds the curvature of T{t), prevents it from 
shrinking to a point, from touching the outer boundary 917, and from undergoing 
topological changes. 

With this property, combining the semiflow for ()1.6|) with the Lyapunov func- 
tional and compactness we obtain the following result. 

Theorem 5.4. Let p > n + 2. Suppose 7 > on (0,itc) and the assumptions 
of (|1.5I) hold true. Suppose that (u, F) is a solution of il.6]) in the state manifold 
SA4 on its maximal time interval [0,^*)- Assume the following on [0,t*); there is 
a constant M > such that 

(i) \u{t)\^2-2,, + |r(i)|^3-3/p <M< 00; 

(ii) < 1/M < u{t) <Uc~ l/M; 

(iii) T{t) satisfies a uniform ball condition. 

Then — 00, i.e. the solution exists globally, and it converges in SM to some 
equilibrium (uoo,roc) G £■ On the contrary, if (u{t),T{t)) is a global solution in 
SA4 which converges to an equilibrium (uH.,r,) in SM as t 00, then properties 
{i)-{iii) are valid. 

Proof. Assume that assertions (i)-(iii) are valid. Then r([0,t,)) C Wp~'^^^ {^l,r) 
is bounded, hence relatively compact in Wp~'^^^' ^{^l,r). Thus we may cover 
this set by finitely many balls with centers Sfe real analytic in such a way that 
dist^y3-3/p-e(r(t),Sj) < 5 for some j = j{t), t € [0,^,). Let = {t & [0,t*) : 
j{t) = k}. Using for each k a Hanzawa-transformation S^, we see that the pull 
backs {u{t, •) o Sfc : i € J^} are bounded in Wp ^^^(f2 \ Efe), hence relatively com- 
pact in Wp "^^^ ^ {^ \ Sfc). Employing now the results in subsection 5.1 we obtain 
solutions (w^,r^) with initial configurations (u(t),V{t)) in the state manifold on a 
common time interval, say (0,t], and by uniqueness we have 

{u\T),T\T))^{u{t + T),T{t + T)). 

Continuous dependence implies then relative compactness of (w(-),r(-)) in SM. 
In particular, = 00 and the orbit (M,r)(R_(_) C SM is relatively compact. 
The negative total entropy is a strict Lyapunov functional, hence the limit set 
F) C SM of a solution is contained in the set 8 of equilibria. By compactness 
a;(u, F) C SM is non-empty, hence the solution comes close to 8, and stays there. 
Then we may apply the convergence result Theorem 15.11 The converse is proved 
by a compactness argument. □ 

6. The Semiflow without Kinetic Undercooling 

In this section we assume throughout 7(5) = for all s > 0, i.e kinetic under- 
cooling is absent. In this case we may apply the results in [T7] and [TU] too, but we 
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have to work harder to apply them. At first we prove (jl.lip as follows. According 
to (11.81) we know that 



Tr{ur)V := (cjr(wr) - n'(r))V = ^^^{divr(drK)Vr?/r) + ld{u)d,uj}. 

Kr(ur) 

Next we observe 
X'iur) 



Kr(ur) 



divr((ir(wr)Vrur) 

- ^ -divr(dr(ur)VrA(wr)) - ^^^A"(7ir)|Vrwrp 



divp(^VrAK)) - ^{a"(.p) - A'(.r)^- 



AF/.r(ur)^|A"(.r) - A'(.r)^}|Vr.r| 



where /ir denotes the antiderivative of drA'/Kr with ft.r(Mm) = 0. We note that 
by a partial integration 

h^{s) = A(.)^ - r A(r)(^)'(r)dr A(s)Mf) _ 
Now employing A(ur) = — ■H(r) leads to the identity 
TT{ur){V-^^n{T)~fr{ur)} 

^ ^Adiu)d.ul - ^{A"(.r) - X'iur f^AWrurl^ 
+ [wr(ur) ~ trLp]ft.r(ur), 
hence applying the inverse of Tr{ur) we arrive at 

Kr{ur)V ~ driur)H{T) = KrK){/rK) + Fr{u, ur)}, (6.1) 

where 

Friu,ur) = [Kr{ur)Triur)]-^X'{ur)ld{u)d,ul 

- drK)[(A"(ur) - X'{ur)K'riur)/Kriur)]\Vrur\^ 
+ Kr(Mr)[wr(ur) - trLp]/ir(Mr)}- 

In the sequel we will replace the Gibbs-Thomson law by the dynamic equation 
?T|) plus the compatibility condition (p(ura) + o'{uto)'H(To) = at time t = 0. 
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6.1. Local Well-Posedness. To prove local well-posedness we employ the direct 
mapping method as introduced in Section 3. As base space we use as in Section 5 

and we let Xi, Xj and X-y,^ as defined there. 

We rewrite system (13. 2p abstractly as the quasilinear problem in Xq 

z + Aoiz)z = F„{z), z(0) = zo, (6.2) 

where z = {v,vr,p) and Zq = (wq, "J^roj Po)- Here the quasilinear part Aq{z) is the 
diagonal matrix operator defined by 

■{d{v)/K{v)){A-M2{p):S/') 



'Ao{z) = diag 



{drivr)/Krivr))iPr{p)Mo{p)f : V| 



idrivT)/KrivT))gip) : V 



with M2(p) = Mi{p) + Mj{p) - Mi{p)Mj{p). The semilinear part Fo{z) is given 

by 

np)v + -^{d'iv)\il - Afi(p))Vz;|' - d{v)({I - M,{p)) : VMi(p)|V«)} 

^{6(i;r, P)v - [{l{vr) + lrivr)H{p) - 7(«r)/3(p)atp]/?(p)9tp + mg} 

{dr{vr)/Kr{vr))T{p) + {fr{vr) + Fr{v,vr, p)} / P{p) 

where by abuse of notation here means the transformed Fj^ introduced previ- 
ously, and where 

TO3 = -drivr)iPr{p)MQ{p))^ : V|t;r -C{vr,p)vr. 

Again, the first two components of Fo{z) contain the time derivative dtp. We 
replace it by the transformed version of (|6.ip 

dtp = {Mvr) + Friv, vr,p) + dr(«r)/«r(t'r)H(p)}//3(p), 

to see that it leads to a lower order term, as in Section 5. 

Provided that Tro(fro) is invertible we may proceed as in Section 5, applying 
Theorem 2.1 in [10 , to obtain local well-posedness, i.e. a unique local solution 

z e «); ^o) n ^^,^((0' «); ^i) ^ Ci[0, a];X^^,,) n C((0, a];X^) 

which depends continuously on the initial value zq G B. The resulting solution 
map [zq I—)- z{t)] defines a local semifiow in X^^^. 

6.2. Nonlinear Stability of Equilibria. Let e, = {ut,ur*,T^,) denote an equi- 
librium as in Section 4. In this case we choose S = F^^ as a reference manifold, and 
as shown in the previous subsection we obtain the abstract quasilinear parabolic 
problem 

i + Ao(z)z = Fo(z), z(0) = zo, (6.3) 
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with Xq, Xi, X^ as above. We set — (w*, ur*, 0). Assuming well-posedness and 
C* 7^ 1 in stability condition, we have shown in Section 4 that the equihbrium e* is 
normally hyperbolic. Therefore we may apply once more |17j . Theorems 2.1 and 
6.1 to obtain the following result. 

Theorem 6.1. Letp > n + 2. Suppose 7 = 0, cr G C^(0,itc), and the assumptions 
of p.Sp hold true. As above £ denotes the set 0/ equilibria of (|5.3p . and we fix 
some e £. Assume that the well-posedness condition 

U ^ and uJl/a^ ^ Kr*{n ~ 1)/RI (6.4) 

is satisfied. Then we have 

(a) //r* is connected and < 1, or if Kr*{ri — l)/i?2 > u*l1/a<, then z* is stable 
in X^, and there exists <5 > such that the unique solution z{t) of h5.S\) with 
initial value zq € Xj satisfying \zo — z*!^ < S exists on R+ and converges at an 
exponential rate in X^ to some z^o & £■ as t ^ 00. 

(b) If KY*{n — 1)/Rl < uJl/a^,, and ifV^ is disconnected or if (^^ > 1 then z* 
is unstable in X-y and even in Xq. For each sufficiently small p > there is 
6 G (0,p] such that the solution z{t) of (j5.3p with initial value zq G Xy subject to 
|zo — z*|^ < 6 either satisfies 

(i) distx.y{z(to); £) > p for some finite time to > 0; or 

(ii) z{t) exists on R+ and converges at exponential rate in X-y to some Zoo G £. 

Thus the only cases which are excluded are = 1, and the two values where 
the well-posedness condition (|6.4p is violated. 

6.3. The Local Semiflow on the State Manifold. We define the state mani- 
folds SMo for ()1.6p in case 7 = as follows. 

SMo ■■= {(M,r) G c{n) X Mn^ -. u g w^-'^'p{n\v), r g w^-^'p, 

< u < Uc in fj, duU = on dO,, (6.5) 
A(-ur) + 'H(r) = on T, Tt{ut) is invertible in L2(r)}. 

Charts for this manifold are obtained by the charts induced by M'H^{Vl) followed 
by a Hanzawa transformation as in Section 3. 

Applying the result of subsection 6.1 and re-parameterizing the interface re- 
peatedly, we see that (|1.6p with 7 = yields a local semiflow on SM.Q. 

Theorem 6.2. Letp > n + 2. Suppose 7 = 0, cr G C^(0,Wc), and the assumptions 
of (jl.Sp hold true. 

Then problem U.6]) generates a local semiflow on the state manifold S Mq . Each 
solution (u,r) exists on a maximal time interval [0,t*), where t,, — t*(uo,ro). 
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6.4. Global Existence and Convergence. In addition to the obstructions to 
global existence for the Stefan problem with variable surface tension in the presence 
of kinetic undercooling there is an additional possibility for loss of well-posedeness: 

• regularity: the norms of u{t) or T{t) become unbounded; 

• well-posedness: the temperature may reach or it j,; or 
Triur) may become non-invertible; 

• geometry: the topology of the interface changes; 
or the interface touches the boundary of fl; 

or the interface contracts to a point. 

We set £o — SAio n £. As in Section 5, combining the semiflow for ()1.6p with the 
Lyapunov functional and compactness we obtain the following result. 

Theorem 6.3. Let p > n + 2. Suppose 7 = 0, cr e C^(0,Wc), and the assumptions 
of (|1.5I) hold true. Suppose that (u, F) is a solution of lll.6\) in the state manifold 
SMo on its maximal time interval [0,i*). Assume the following on [0,^*); there 
is a constant M > such that 

(i) \u{t)\^,-.,„ + |r(i)|^3-3/p < M < cx); 

(ii) < 1/M < u{t)) <Uc~ 1/M; 

(iii) > 1/AI holds for the eigenvalues of TY{t){ur); 

(iv) T{t) satisfies a uniform hall condition. 

Then t^ = 00, i.e. the solution exists globally, and it converges in SMq to an 
equilibrium (uoo,roo) G £0 ■ Conversely, if (u(t),T{t)) is a global solution in SMq 
which converges to an equilibrium {uca,^oo) G £q in SMq as t ^ 00, then the 
properties {i)-{iv) are valid. 

Proof. The proof follows the same lines as that of Theorem 15.41 □ 
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